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^ Over the last number of years spectroscopic studies have strongly supported 

^ the assertion that protostellar accretion and outflow activity persists to the lowest 

masses. Indeed, previous to this work the existence of three brown dwarf (BD) 
^ outflows had been confirmed by us. In this paper we present the results of our 

l/-^ latest investigation of BD outfiow activity and report on the discovery of two 

^ new outfiows. Observations to date have concentrated on studying the forbidden 

CN emission line (FEL) regions of young BDs and in all cases data has been collected 

O using the UV- Visual Echelle Spectrometer (UVES) on the ESO VLT. Offsets 

in the FEL regions are recovered using spectro-astrometry. Here ISO-Oph 32 
O is shown to drive a blue-shifted outflow with a radial velocity of 10-20 kms~^ 

_ ^ and spectro-astrometric analysis constrains the position angle of this outflow to 

^X^ 240° ± 7°. The BD candidate ISO-Chal 217 is found to have a bipolar outflow 

bright in several key forbidden lines {Y rad = -20 kms^^, +40kms^^) and with 

a PA of 190°-210°. A striking feature of the ISO-Chal 217 outflow is the strong 
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asymmetry between the red and blue-shifted lobes. This asymmetry is revealed in 
the relative brightness of the two lobes (red-shifted lobe is brighter), the factor of 
two difference in radial velocity (the red-shifted lobe is faster) and the difference 
in the electron density (again higher in the red lobe). Such asymmetries are 
common in jets from low mass protostars (0.1 Moto 2 M0) and the observation 
of a marked asymmetry at such a low mass ( < 0.1 Mq) supports the idea that BD 
outflow activity is scaled down from low mass protostellar activity. Also note that 
although asymmetries are unexceptional, it is uncommon for the red-shifted lobe 
to be the brightest as some obscuration by the accretion disk is assumed. This 
phenomenon has only been observed in one other source, the classical T Tauri 
(CTTS) star RW Aur. The physical mechanism responsible for the brightening 
of the red-shifted lobe although as yet unknown must also now apply at BD 
masses to include the ISO-Chal 217 outflow. In addition to presenting these 
new results, a comprehensive comparison is made between BD outflow activity 
and jets launched by CTTSs. In particular, the application of current methods 
for investigating the excitation conditions and mass loss rates (Mout) in CTT jets 
to BD spectra is explored. Where possible the Bacciotti & Eisloffel technique is 
used to study the ionisation fraction, electron temperature and total density. For 
LS-RCrA 1, ISO-Chal 217 and ISO-Oph 102 Mout is measured to be in the range 
10~^° to 10^^ M0yr~^ using a method based on the luminosity of the [OI]A6300 
and [SII]A6731 lines. Mass loss rates for our sample of BD outflows are found to 
be comparable to the mass accretion rates. Overall, as our results and discussion 
show, what is currently known about outflow activity in the BD mass regime 
points to strong similarities between BD outflows and CTT jets. 

Subject headings: — stars: low mass, brown dwarfs — stars: formation — ISM: 
jets and outflows 

1. Introduction 



The outflow phenomenon, much studied in young stars has been observed in a wide 
variety of astrophysical objects from low mass protostars to active galactic nuclei (AGN) 



(Reipurth & Bally 


2001 


Camenzinc 


I 200^ 


)). Brown Dwarfs (BDs) are the newest objects to 


be added to this list ( 


Fernandez & Comeron 


2001 


Barrado y Navascues et al. 


2004 


Whelan 


et al. 


2005 


2007 


Phan-Bao et al. ! 


2008 


Whelan et al. 


2009 


). We are currently leading a 



project to search for and study outflows driven by young BDs and very low mass stars. 
The overall aim of this work is to explore outflow activity in the BD mass regime through 
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comparison with observations of outflows driven by low mass protostars and the classical T 
Tauri stars (CTTSs) in particular. Much emphasis has been placed on understanding how 



accretion activity (Gatti et al. 2006 Herczeg et al. 2009) and planet formation processes 



(Pascucci et al. 2009) differ between low mass protostars and young BDs. However, little is 



still known about whether the outflow mechanism varies as mass decreases towards the BD 
regime. Observations aimed at detecting outflows from such low mass objects are difficult 
and rare and to date six candidates have been investigated by us spectro-scopically at visible 
wavelengths. 

Outflows from low mass protostars are revealed in several forms, from slow wide angled 



winds to molecular outflows and jets (Reipurth & Bally 2001). A protostellar jet is a rel- 
atively fast (typical radial velocity = 200 kms~^) collimated outflow normally seen within 
a few 1000 AU of a young star and bright in optical and/or near-infrared shock tracers 
such as [SII]A6731, [Fell] at 1.644/xm or H2 at 2.122 /xm. CTTSs are Class II low mass 
protostars ( Lada|[l987 Andre et al. 1993), meaning that they are coming towards the end 



of their accretion/outflow phase and have crossed the so-called "birth-line" (Stabler 1983). 
As they are optically visible the jet can be traced very close to the protostar, hence, they 
offer a unique opportunity to study the launching of protostellar jets at high spatial reso- 



lution (Ray et al. 2007). Jets from CTTSs are traditionally probed at forbidden emission 



line (FEL) wavelengths and long-slit and integral field spectroscopic techniques have been 



hugely important in their study (Dougados et al. 2000; Whelan et al. 2004). That CTT-like 



outflows could be launched by actively accreting BDs was first suggested when high quality 
spectra revealed the presence of FEL regions in the optical spectra of BDs known to be 



accretors (Fernandez & Comeron 2001). While the FEL regions were easily identified they 



were considerably fainter than those detected in the spectra of CTTSs (Masciadri & Raga 



2004 ) and any extension in the form of an outflow was not apparent ( Fernandez & Comeron 



2001, 2005). Hence their origin in an outflow could not be established. The critical density 
region for FELs occurs much closer to the driving source of the outflow for nearby (~ 150 
pc) young BDs than for low mass protostars and for BDs is estimated to be on scales of ~ 



100 mas (Whelan et al. 2005). As the most intense forbidden emission coincides with the 



critical density region it is currently challenging to directly resolve a BD outflow in FELs 
and this goal has not yet been achieved. Our approach to this problem has been to obtain 
high quality spectra using the UV- Visual Echelle Spectrometer (UVES) on the European 
Southern Observatory's (ESO) Very Large Telescope (VLT) and to recover the spatial offset 
in the region of forbidden emission using spectro-astrometry. 

The sources targeted so far meet at least one of two criteria. Firstly they all have masses 
at or significantly below the hydrogen burning mass limit (HBML) and secondly they are 
known to be powerful accretors. Their classification as a strong accretor is based on studies 
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of their Ha lines ( Jayawardhana et al. 2003a) and on measurements of their mass accretion 
rates (see Table 6). In the vast majority of cases the mass of young BD candidates cannot be 
directly measured and estimates are based on theoretical predictions (Luhman et al. 2000) 
and/or spectral modelling (Mohanty et al. 2004). The objects investigated by us so far are 
listed in Table[T]along with their theoretical masses and spectral types. Previous to this paper 
we have reported on the discovery of and discussed the optical outflows from, ISO-Oph 102 
dWhelanet aL]|2005| , 2MASS1207-3932 ( |Whelan et aL][20^ and LS-RCrA 1 ( |Whelan et al 



2009). With a mass of only 24 Mjup 2MASS1207-3932 is at present the lowest mass galactic 



object driving an outflow. In this paper we present the results of our recent observations 
of ISO-Chal 217, DENIS-P J160603.9-205644 and ISO-Oph 32 and report for the first time 
outflows driven by ISO-Chal 217 and ISO-Oph 32. Of the six objects observed by us to date 
only one, ISO-Chal 217 has a theoretical mass which is just above the HBML, placing it 
at present just outside the mass range for a BD. However, with a predicted mass of just 80 
^jup it is still one of the lowest mass objects known to launch an outflow. Therefore for the 
purpose of the discussion we will include it in this paper with the likely BDs. In addition to 
presenting our spectro-astrometric analysis, a technique for constraining the position angle 
(PA) of the BD outflows is outlined, methods for studying the physical conditions of the BD 
outflow and of estimating the mass outflow rate are discussed and similarities between BD 
outflows and CTT jets highlighted and considered. 



2. Observations and Spectro-astrometric Analysis 



Spectra (at orthogonal slit positions) of ISO-Chal 217, ISO-Oph 32 and DENIS-P 
J160603. 9-205644 were obtained with the high resolution spectrometer UVES (CD3 dis- 
perser, spectral range 4810-6740 A, R=40,000;), on the ESO VLT UT2 ( [Dekker et aL]|2000[ ), 
as part of the observing program 079.C-0375(A). In each case the exposure time was 2700s, 
the slit width 1" and the seeing varied between O'.'S and 0'.'9. Spectra were analysed us- 



ing spectro-astrometry as outlined in previous papers (Whelan et al. 2005 2007). Also see 



Whelan & Garcia (2008) for details on basic data reduction steps and continuum/night sky 
line subtraction. To summarise spectro-astrometry involves the measurement of the spatial 
centroid of an emission line (using Gaussian fitting) as a function of velocity. The centroid 
of the continuum emission is taken as the position of the star/BD and the emission line 
centroid is measured with respect to this position. 

Problems specific to the application of spectro-astrometry to BD spectra are the lack of 
or weakness of continuum emission from the BD in the range 5000-8000 A and the faintness 
of the line emission itself. In order to overcome these problems and increase the spectro- 
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astrometric accuracy, the spectra are summed or smoothed in the dispersion direction. In 

1/2 

the absence of detector noise, the spectro-astrometric accuracy is dependent on Np' where 
Np is the number of detected photons in the hne/continuum (Whelan et aL 2005). In the 
cases where spectra are summed, the hne and continuum are summed separately and by 
differing amounts, so that the spectro-astrometric error in the line and continuum position 
is comparable. Recall that the continuum must be removed before the line centroid is fitted. 
In these cases, points are a moving sum across the line and continuum (see Figures 3-7 
and Figure 11). For BDs where the line emission is weaker than the underlying continuum 
emission it was more effective to decrease the spectral resolution and hence increase the 
signal by smoothing the spectrum using a Gaussian filter ( |Whelan et aL]|2007D . A gain this is 
done so that errors are comparable. The [OI]A6300 line region in the spectra of ISO-Oph 32 
was smoothed using a Gaussian filters of fwhm = 0.2 A x 0'.'32. The results of the spectro- 
astrometric analysis of the smoothed spectra are shown in Figures 9 and 10. In all cases the 
dashed line delinates the 1-a error in the continuum (and hence line) centroid. Finally, the 
issue of spectro-astrometric artifacts must be addressed. Spectro-astrometric artifacts, or 
false spectro-astrometric signals, can be introduced into a spectrum for a variety of reasons 
(Podio et al. 2008). In our studies of BD outflows we have ruled out the influence of such 
artifacts by demonstrating that lines like HeIA6678, which will not form in an outflow, do 
not have any spectro-astrometric signal. If the outflow signatures detected in the FELs were 
due to artifacts then these lines would also show similar offsets. In addition, note that the 
Ha line in the majority of cases does not show any offset. This is in agreement with the 
idea that most of the Ha emission comes from the accretion zone close to the source. For 
all the BDs, spectra are presented at several PAs, and offsets measured at the different PAs 
agree. Lastly the position velocity (PV) diagram of the emission line regions of ISO-Chal 
217 supports the spectro-astrometric analysis. 



3. Results of Spectro-astrometric Analysis 



In this section the results of the spectro-astrometric analysis of the FE and Ha line 
regions of, ISO-Chal 217, ISO-Oph 32 and DENIS-P J160603.9-205644, are presented. The 
particulars of the spectro-astrometric technique, as applied to BD spectra, are discussed 
above. As well as recovering the spatial offset in the FEL regions of our BD candidates 
(FELs were only detected in the case of ISO-Chal 217 and ISO-Oph 32) this technique is 
also used to estimate the PA of the new BD outflows, through comparison of the spectro- 
astrometric offsets measured at orthogonal slit positions. Jets driven by CTTSs were first 



observed using narrow band imagery and long slit spectroscopy (Mundt & Fried 1983; Solf 



fc Boehm|[l993| . |Hirth et al.| ( |I994| [l997| estimated the micro-jet PA for a sample of CTTSs 
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by taking spectra at several slit PAs and plotting the positional offset and spatial full width 
half maximum (FWHM) of the [SII]A6731 line as a function of slit PA. They argue that the 
[SII]A6731 line is a better indicator of the outflow direction than the [OI]A6300 line, as its 
lower critical density means that it traces the jet at a greater distance from the source. Closer 
to the source other outflow components e.g the low velocity component could influence any 



measurement of the jet PA. Although, the approach of Hirth et al. (1994) would provide a 
more straight-forward way of estimating the outflow PA for our BD candidates, the number 
of spectra that could be collected was limited by the faintness of the FELs and thus the 
on source time needed to detect them with a reasonable S/N (on average 45 mins with the 
VLT). For UVES, a negative spectro-astrometric offset, at a N-S slit PA equals a southerly 
displacement in an emission feature and a positive offset a northerly displacement. For an 
E-W slit P.A, negative and positive offsets correspond to westerly and easterly directions 
respectively. This is illustrated in Figure 1. As demonstrated in Figure 2. comparison 
of the signs of the N-S, E-W spectro-astrometric offsets allows the quadrant in which the 
blue-shifted outflow lies to be ascertained. Once this is known the actual outflow PA (E of 
N) can be derived. The dashed, arrowed lines represent the blue-shifted lobe of the outflow 
as located in each quadrant. The angles Ai to A4 are simply estimated from tanfAi ^) = 
(offset for E-W sUt)/ (offset for N-S sUt). 

Error is introduced by the fact that the positional offset in the FEL regions is not directly 
resolved but recovered using spectro-astrometry. The accuracy which one can measure the 
positional displacement of an unresolved emission feature using Gaussian fitting depends 
primarily on the S/N of the emission feature. However the displacement is also weighted by 
the uniformity of the emission. One emission feature could actually be made up of several 
unresolved features of differing intensity. Where the spectro-astrometric technique measures 
the spatial mean of the feature to be, is influenced by where the maximum intensity feature 
lies. In other words an unresolved jet will be measured to be less extended than it really is 
(using spectro-astrometry), if it actually consists of bright inner knots and lower brightness 
material further out. For the above reasons the estimates of the PAs of the ISO-Chal 217 and 
ISO-Oph 32 outflows given here are a first approximation and the errors quoted are taken 
from the spectro-astrometric uncertainity (dependent on S/N). Also note that outflow PAs 
measured here only apply to the outflow as studied close to its driving source. It is possible 
that some bending may occur in the outflow at large distances which will change the outflow 



PA (Fendt &; Zinnecker 1998). However we expect any such variation to be within the 
errors quoted in Table 3. Below, the results for ISO-Chal 217, ISO-Oph 32 and DENIS-P 
J160603. 9-20564 are presented separately. Tables [2] and [3] give the FELs detected for each 
source along with the radial velocity (measured with respect to the systemic velocity), offset 
measured and estimate of the outflow PA. 
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Source 


RA (J2000) 


Dec (J2000) 


Spectral Type 


Mass (Mjup) 


ISO-Chal 217 


11 09 52.0 


-76 39 12.0 


M6.2 


80^ 


2MASS1207-3932 


12 07 33.4 


-39 32 54.0 


M8 


242 


DENIS-P J160603.9-205644 


16 06 03.90 


-20 56 44.6 


M7.5 


40^ 


ISO-Oph 32 


16 26 22.05 


-24 44 37.5 


M8 


40* 


ISO-Oph 102 


16 27 06.58 


-24 41 47.9 


M6 


60* 


LS-RCrA 1 


19 01 33.7 


-37 00 30.0 


M6.5 


35-72^ 



Table 1: The spectral type and predicted mass of the BD candidates investigated by us 
to date. All sources except DENIS-P J160603. 9-205644 are found to drive outflows. The 
numbers 1-6 refer to the papers giving the estimates of the mass of each source, where 
1 



MuzeroUe et al. (2005), 2=Mohanty et al. (2007), 3=Mohanty et al. (2004), 4=Natta et 



al. (2002) and 5=Barrado y Navascues et al. (2004) 



Object 


Slit PA (°) 


Line 


V (kms-i) 


Offset (mas) 


ISO-Chal 217 





[OI]A6300 


-20/30 


-50/50 (± 6) 






[OI]A6363 


-10/28 


-50/45 (± 9) 






[SII]A6716 


-20/42 


-180/200 (± 33.5) 






[SII]A6731 


-20/38 


-180/180 (± 33.5) 




90 


[OI]A6300 


-10/20 


-30/25 (± 6) 






[OI]A6363 


-10/20 


-40/10 (± 12) 






[SII]A6716 


-15/40 


-46/40 (± 14) 






[SII]A6731 


-20/45 


-75/55 (± 14) 


ISO-Oph 32 





[OI]A6300 


-10 


-60 (± 10) 




90 


[OI]A6300 


-20 


-105 (± 19) 



Table 2: The FELs of ISO-Chal 217 and ISO-Oph 32 detected at slit PAs of 0° and 90°. 
Also given are the radial velocity of the lines (blue and red-shifted components in the case 
of ISO-Chal 217) and the offsets in milli-arcseconds recovered using spectro-astrometry. An 
estimate of the outflow PA is made by comparing the offsets measured at the orthogonal 
PAs. 
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N 



UVES slit at P.A. = degrees 



it 



+ve offset, emission extended in northerly direction 



i 



wavelength 



UVES slit at P.A. = 90 degrees 



★ 



+ve offset, emission extended in easterly direction 



i 



wavelength 



-ve offset, emission extended in soutlierly direction 



-ve offset, emission extended in westerly direction 



Fig. 1. — Figure illustrating the direction on the sky which the sign of the spectro-astrometric 
offsets correspond to, for different slit PAs. Above left, the slit is positioned at a PA of 
0°. Below left, the resulting spectrum is shown as a strong emission line feature (assumed 
to originate in an outflow) superimposed on the continuum emission. A positive spectro- 
astrometric offset corresponds to extended emission in a northerly direction and a negative 
offset a southerly direction. In the case of an E-W slit position (PA=90°) a positive offset 
corresponds to an easterly direction and negative to a westerly. Wavelength increases towards 
the left. 
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+ve offset N-S slit P.A. 
+ve offset E-W slit P.A. 

Outflow P.A. (E of N) = A2 



-ve offset N-S slit PA. 
+ve offset E-W slit PA. 

Outflow RA. (E of N) = I8O-A3 



N 



A2 



Ai 




A3 



A4 



^ +ve offset N-S slit PA. 
..• -ve offset E-W slit RA. 

Outflow RA. (E of N) = 180-Ai 



-ve offset N-S slit PA. 
-ve offset E-W slit PA. 

Outflow RA. (E of N) = I8O+A4 



Fig. 2. — Figure illustrating how the signs of the spectro-astrometric offsets can be used to 
determine the quadrant in which the outflow lies. The dashed lines represent the blue-shifted 
outflow. Angles A1....A4 are estimated from the following formula tan(A) = (offset for E-W 
slit) /(offset for N-S slit). As negative offsets are measured for both the ISO-Chal 217 and 
ISO-Oph 32 outflows the blue-shifted lobes lie in the fourth quadrant. 



3.1. ISO-Chal 217 



ISO-Chal 217 was first reported as part of an ISOCAM survey of tlie Cliamaeleon I 



dark cloud (Persi et al. 2000). Current estimates place its mass at 80 Mjup (MuzeroUe et 



aL] |2005D and its spectral type between M7.5 and M6.25 ( [Lopez Marti et al.||2004| [MuzeroUe 



et al. 2005). Hence, it currently lies on the boundary between being a BD or an ultra- low 



mass star. The accretion properties of ISO-Chal 217 have been much studied in recent years. 



Apai et al. 


(2005 


) confirm mi 


and 


Scholz & Jayawardhana 



significant variability in its accretion activity. Scholz & Jayawardhana (2006) also observe 
signatures of outflow activity in their spectra of ISO-Chal 217. These signatures are in 
the form of FELs of S. Here the origin of these (and other) FELs in an outflow driven by 
ISO-Chal 217 is confirmed for the first time. 

UVES spectra of ISO-Chal 217, taken at orthogonal slit positions were obtained on 
May 1st 2007. Strong [SII]AA6716,6731, [OI]AA6300,6363 and Ha emission was detected 
at both PAs along with weak [NII]A6548 emission. Blue and red-shifted components to the 
[SII]AA6716,6731 and [OI]AA6300,6363 lines were detected at an average radial velocity (w.r.t 
the systemic velocity) of -16 and 32 kms^^. The systemic velocity of ISO-Chal 217 is taken 



James et al. 


2006; 


Bary et al. 


2008 



and Ha lines were analysed with spectro-astrometry and the analysis revealed the blue and 
red-shifted emission to be offset in opposing directions. Hence, ISO-Cha I 217 is driving a 
bipolar outflow (see Figures 3, 4, 5, 6, 7). As is the case in CTTS the [SII]AA6716,6731 lines 
trace the outflow further from the source than the [OI]AA6300,6363 lines (refer to Table [2]). 
Emission is clearly more extended at 0° than at 90° and a comparison between the two is 
used to constrain the outflow PA to a range of 193° to 208° (E of N). Note that, the lines 
which trace the outflow closest to the source give a shghtly larger PA on average than the 



Object 


Line 


Outflow PA (°) 


ISO-Chal 217 


[OI]A6300 


209 (± 8) 




[OI]A6363 


206 (± 15) 




[SII]A6716 


193 (± 5) 




[SII]A6731 


200 (± 6) 


ISO-Oph 32 


[OI]A6300 


240 (± 7) 



Table 3: Table giving the position angles of the two outflows observed. The PAs are estimated 
from orthogonal spectra as described in Section 3. For ISO-Chal 217 a bipolar outflow is 
uncovered and estimates given above are an average of measurements made for the red and 
blue-shifted lobes. For ISO-Oph 32 only the blue-shifted outflow is detected. 



more displaced line regions (Table 3), although all values agree within the errors. Analysis 
of the Ha line reveals no detectable contribution from the outflow. Clearly the Ha line is 
coincident with the continuum position at both position angles. 

A significant feature is the obvious asymmetry in the outflow as traced by [SII]. The 
red-shifted lobe is much brighter in [SII] AA6716, 6731 than the blue-shifted and there is a 
velocity difference of a factor of 2 between the red and the blue lobe. This is especially 
apparent in the case of the brighter [SII]A6731 line at a PA of 0° (Figure 7). Although 
velocity and intensity asymmetries are commonly seen in CTT jets (Hirth et al. 1994) and 
while several case have been documented where the red-shifted lobe is the faster it is unusual 
for the red-shifted lobe to be the brightest. In the majority of CTT jets the circumstellar disk 
obscures the red-shifted flow. The CTTS RW Aur offers the most striking example of such 
an asymmetry (Hirth et al. 1994 Lopez-Martm et al. 2003). In this case a similar velocity 
asymmetry has been reported but the slower lobe (red) is the brightest. The ISO-Chal 217 
asymmetry is clearly visible in the PV diagrams presented in Figure 8. For the purpose 
of the PV diagrams the [SII]A6731 line region was smoothed, using an elliptical Gaussian 
filter. The PV diagrams support the spectro-astrometric analysis as the opposing offset in 
the line region at 0°, measured with spectro-astrometry to be ~ 180 mas, is resolved in the 
corresponding PV plot. In the PV plot at 90° the line appears coincident with the source 
position again consistent with the much smaller offsets measured at this slit position (~ 50 
mas). 



3.2. DENIS-P J160603.9-205644 



DENIS-P J160603. 9-205644 was identified as a BD as part of a survey of the Upper 



Scorpius OB association Martin et al. (2004). Mohanty et al. (2005) classify it as an accretor 



(Ha 10% width is given as 306 kms ^) and give its spectral type and mass as M7.5 and 40 
Mjup. However, no FELs or additonal signatures of outflow activity are found in the spectra 
of DENIS-P J160603.9-205644. The spectro-astrometric analysis of the Ha line at 0° and 90° 
is shown in Figure 11. While the line profile is asymmetric the emission is coincident with 
the continuum and therefore must originate at the star. The lack of a spectro-astrometric 
signal in the Ha line rules out any complications with spectro-astrometric artifacts. Although 
DENIS-P J160603. 9-205644 is classified as an accretor it is likely that it is the least accreting 



object in our sample. Mohanty et al. (|2005|) give upper limits to its mass accretion rate based 

'acc for 



Herczeg et al. (2009) measure 



on the non detection of the Ca II 8662 A feature 
this object at 2.5 x 10~^^ M0yr~^ consistent with the upper limits set by Mohanty et al. 



(2005). Also note that Herczeg et al. (2009) do detect some faint forbidden emission in the 
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spectrum of DENIS-P J160603. 9-205644. Their observations were made with LIRIS on Keck 
I at a spectral resolution of R=1400. Given the relatively high spectral resolution of our 
observations and the small accretion rate of this object it is not surprising that no detection 
is made here. 



3.3. ISO-Oph 32 



Natta et al. (2002) observed this source in the near-infrared and estimated its spectral 
type and mass at M7.5 and 30-50 Mj^/p. They also argue that the mid-infrared excess 
detected in this object is consistent with the presence of an irradiated disk. In a later study 



Natta et al. (2004) added optical spectra to their near-infrared data. In particular they use 



the Ha line as a test of accretion activity. In the both the 0° and 90°, high resolution UVES 
spectrum presented here the [OI]A6300 line is detected at at S/N of just above 4. As the 
[OI]A6363 line is 3.5 times fainter than the [OI]A6300 line, it is marginally detected at both 
PAs and therefore not analysed with spectro-astrometry. The lower critical density lines are 
below the detection limit in our spectra. As shown in Figure 9 analysis reveals the [OI]A6300 
line to be offset in a southerly direction to ~ 60 mas at 0° and ~ 105 mas at 90°. The average 
radial velocity (w.r.t the systemic velocity) of the two lines is 15 kms~^ and comparison of 
the measured offsets at the orthogonal slit positions constrains the PA of the blue-shifted 
outflow at ~ 240°. The systemic velocity of ISO-Oph 32 is taken as +7 kms~^ by Whelan 



et al. (2005). Again no spectro-astrometric signal is detected in the Ha line, confirming the 



origin of the bulk of the emission in the accretion flow. 



4. Measuring the Mass Outflow Rate for BD Outflows 



at 



Numerous observational studies have constrained the mass outflow rate [Mout] in CTTSs 
1-10 % of the mass accretion rate ( Hartigan et al.|[l994 Podio et al. |2006 Garcia Lopez 



et al. 



2008). The ratio between the mass outflow and accretion rates {Mout/ Mace) in young 



stars is a deflning characteristic of the mechanism which is at work to launch protostellar 



outflows (Ferreira et al. 2006 Ray et al. 2007). If a similar ratio between outflow and infall 



is measured for accreting BDs it will provide compelling evidence that the outflows we are 
beginning to detect at very low and sub-stellar masses are in fact scaled-down from low mass 
protostellar outflows. The spectra of low mass protostars contain a wealth of forbidden and 
permitted lines whose excitation properties are well known. Hence, the ratios and intensities 
of these lines can be used to probe the physical conditions of the emitting regions. Spectral 
analysis of known outflow tracers, typically in the optical and near-infrared regimes, has 
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allowed Mout to be estimated for low mass protostars. In comparison, the spectra of young 
BDs, have been little investigated. It is only in the last number of years that optical FELs 
have been detected in BD spectra and comparatively little work has been done in the near- 
infrared. 

Two distinct methods have been used to measure Mout in low mass protostars. Here the 
application of these methods to the problem of measuring mass outflow rates at BD masses, 
is discussed. The first method (hereafter method A) relies on an estimate of the total density 
in the outflow uh, which in turns requires the electron density (ug) and ionisation fraction 
(xg) to be known ( Nisini et al.||2005 Podio et al.||2006 ). These plasma parameters can be cal- 
culated using the Bacciotti & Eisloffel (BE) technique ( Bacciotti fc Eisloffel||1999a Bacciotti 
et al. 1999b). The BE technique employs selected optical transitions of S+, O" and N+, to 



derive Xg and Tg (electron temperature). The electron density is estimated from the ratio of 
the [SII] lines, [SII]A6716/[SII]A6731 ( [Osterbrock |[I994| and once this quantity is known it is 



combined with Xe to derive uh- Podio et al. (2006 ) adopt the BE technique to study the phys- 



ical conditions in a sample of protostellar outflows and thus calculate their mass loss rate. 
As the BE technique calculates Tg and from the [Nil] AA (6548+6583)/ [01] AA (6300+6363) 
and [SII]AA(6716+6731)/[OI]AA(6300+6363) line ratios, it requires the use of a given set of 



elemental abundances. Podio et al.| (2006) investigate any variations in estimates of Xg and 
Te for four different abundance sets and calculate Mout using. 



Mout = finiHnHiirjVj (1) 

where /i=1.41 is the mean atomic weight, m^^ the mass of a proton, hh the total density 
and Yj, Yj the jet radius and velocity respectively. Here we firstly discuss the analysis of our 
BD spectra with the BE technique. 

The BE technique is easily apphed to the spectrum of a typical CTTS where the required 
FELs are well resolved spatially and spectrally with a high S/N. The line flux can be measured 
using several different line fitting techniques. The main source of error comes from the 
uncertainty in the measurement of the line fluxes which depends on the S/N ratio of the line. 
In order to use the BE technique to investigate a BD outflow the [OI]A6300, [SII]A6716, 6731 
and [NII]A6583 lines must all be detected at a sufficient S/N. The [OI]A6363 and [NII]A6548 
lines can be substituted for the [OI]A6300, [NII]A6583 lines, if for example the brighter 
[OI]A6300 line is not covered by the spectral range of the instrument or the [NII]A6548 is 
blended with the Ha line. Both the [OI]A6300 and [NII]A6548 lines can be exchanged with 
the fainter line of their pairs, as the ratio of their intensities is always three. This is due to 
the fact that they come from the same upper level and that the ratio of the A coefficients for 
spontanoeus emission is three. Of the five BD candidates found so far to have outflows only 
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one, LS-RCrA 1 (outflow discussed in Whelan et al. 2009), meets the requirements of the BE 



technique. For the two lowest mass BDs, ISO-Oph 32 and 2MASS1207-3932 (Whelan et al. 



2007[ ) only the [OI]AA6300,6363 lines are observed. An added difficulty for 2MASS1207-3932 
is that while both the blue and red-shifted lobes of the outflow are detected in the [OI]A6300 
line they are not well resolved spectrally hence the emission from both lobes cannot be 
properly separated. While both the [OI]AA6300,6363 and [SII]A6731 lines are found in the 



spectrum of ISO-Oph 102 (Whelan et al. 2005) the [SII]A6716 line is not detected and 



therefore there is no estimate of Ug. Finally for ISO-Chal 217 , while the [OI]A6300 and 
[SII]AA6716,6731 lines are all observed the [NII]A6548 line is not present and the [NII]A6583 
line is too faint to allow the use of the BE technique. However, using the ratio of the 
[SII]AA6716,6731 lines Ug in the blue and red-shifted lobes of the ISO-Chal 217 outflow is 
measured (see Table lil). 



The results of the application of the BE technique to the blue-shifted outflow driven by 



LS-RCrA 1 are reported in Table |4j As described in Whelan et al. (2009) UVES spectra 
of LS-RCrA 1 were taken at several P.As between 45° and 82° and analysis of this data 
constrained the PA of the LS-RCrA 1 outflow at ~ 15°. Using the BE technique Ue, Xg, 
Te and uh are derived at each PA. The quoted errors come from the uncertainty in the 
measurement of the line fluxes. In the case of low mass protostars, values of Xe in the range 



0.01 - 0.4, Te in the range 8000-20000 K and ylh up to 10^ cm ^ are typically measured (iRay 



et al. 2007). The estimates of Ue for LS-RCrA 1 and ISO-Chal 217 are reasonable when 
compared with estimates made for CTT jets and the large errors (up to ~ 40%) are due 
to the low S/N in the [SII] lines, which is a low as 3 in some cases. [Comeron et al. (2003) 



quoted a value of Ue = 2300 cm'^ for LS-RCrA 1. Note that the ISO-Chal 217 red-shifted 
outflow has a higher electron density that the blue-shifted which is a further indication of the 
asymmetry between the two lobes of the outflow. The low S/N of the LS-RCrA 1 [NII]A6583 
line leads to errors of up to 50% in the derived values of Xg and hh- The estimates of Xg 
are very low compared to what one measures for CTTSs and decrease as the slit PA moves 
away from the outflow PA. In addition, values of Tg and iih are very high. With such high 
densities one would expect more efficient cooling and therefore lower temperatures. The 
peculiarities in Xg, Tg and noted here have been observed in CTT flows where the very 
beginning of the flow (0'.'2 to 1") is being tested with BE technique. In the case of the BDs 



we are probing even closer to the driving source of the outflow. Bacciotti et al. (1999b) argue 
that the physical reason for this behaviour could be related to the event that creates the 
ionisation of the flow e.g. a standing shock. While these preliminary results are interesting, 
it is impossible at present to say if this approach can reveal anything meaningful about the 
ionisation fraction / electron temperature and thus shocked emission in BD outflows. Higher 
S/N spectra of a bigger sample of young BD outflows, aimed at resolving the [Nil] lines in 
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particular, should be the next step. 



The value of uh estimated from the above method is used in Equation 1 to derive Mout- 
As the slit PA of 45° lies closest to the estimated outflow PA (15°) nn derived for this slit 
PA is taken. One of the biggest difficulties with measuring Mout for BD outflows is that the 
outflow width and velocity is not known. Here reasonable estimates of tj, vj for LS-RCrA 



1 are made. Ray et al. (2007) show a plot of jet width as a function of projected distance 



from the source. Assuming a similar relationship for a young BD outflow and considering 
that offsets of up to 150 mas (~ 20 AU) have been measured for the LS-RCrA 1 FELs, 10 
AU is taken as a reasoned approximation of the radius of the LS-RCrA 1 outflow. This 
would suggest an opening angle close to 30° for this outflow and is likely an upper limit on 
the outflow radius. The radial velocities of the LS-RCrA 1 FELs are low (~ 25kms^^) and 
at present the inclination angle of the outflow is poorly constrained. Initial observations 
suggested that the low radial velocity of the FELs and the faintness of the source could be 
explained by an edge-on accretion disk and thus an outflow inclined close to the plane of 
the sky. However the lack of red-shifted emission and the evidence for separation between 



low and high velocity emission in the outflow, presented in Whelan et al. (2009), refute this 



argument. Assuming that the low radial velocity measured in the outflow tracers suggests a 
low outflow inclination angle (20-30°), for the purpose of calculating Mout, vj is taken at 75 
kms~^. Combining ylh with ij and vj yields Mout = 2.4 x 10~^ MQyr"^. 

The second widely used method for measuring Mout (hereafter method B) is based on 
the observed luminosity L{line) of an optically thin line such as [SH] or [01]. The luminosity 
of a speciflc line can be derived as follows (Comeron et al. 2003[ ) 



L{Lq) = 6.71 X 10-^D\pc)EW{A)10-^-^^° 



(2) 



(D is the distance to the object, EW the equivalent width of the line and Rq the dereddened 
magnitude) and used to evaluate the mass of the flow M (in the aperture) according to the 



following equations (Hartigan et al. 1995), 



M = 9.61 X 10" 



M = 1.43 X 10" 



:i + !^)(^)Me 
:i + ^)(^)Me 



Up 



(3) 



(4) 







M, 



out 



MVtan/han whcrc Uc, Vtan and \tan are the critical density, outflow tangential veloc- 
ity and the size of the aperture in the plane of the sky. As Ug for CTTSs is comparable to 



the critcal density of [SII]A6731 Hartigan et al. (1995) remove the critical density depen- 
dence from Equation 4 above. For the BD outflows investigated to date Ue < Uc hence this 
dependence is included here. 
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In order to measure L{line) and apply the above method to our sources, Rq, and Vtan 
must be known. Here we estimate M^ut for LS-RCrA 1 , ISO-Chal 217 and ISO-Oph 102. 
The values of EW, Rq, Ue and Vtan used in each case are given in Table |5} For LS-RCrA 1 
Mont is measured at 6.1 x 10"^° Moyr"^ and 2.0 x 10"^° Moyr"^ from the [OI]A6300 and 



[SII]A6731 lines respectively. Vtan is taken at -43 kms~^. This estimate is based on a radial 
velocity of 25 kms~^ and an outflow inclination angle of 30° (see above). For the ISO-Chal 
217 bipolar outflow we chose to estimate Mout from the [SII]A6731 line luminosity. As the 
outflow is resolved spectrally in the [SII]A6731 line (the lobes are blended spectrally in the 
[OI]A6300 line) it was possible to estimate Mout in both lobes of the outflow. For the faster 
(red-shifted lobe) Mout = 3.1 x 10""^" MqJt~'^. This is almost twice the outflow rate of the 
blue lobe, which is measured at 1.8 x 10^^° MQyr"^. A similar difference in Mout is observed 



in the two lobes of the RW Aur asymmetric jet Woitas et al. (2002). Rq for ISO-Chal 217 is 



taken from the R = 19.59 mag dLopez Marti et al.||2004[ ) and Aj given by iLuhmanl (|2004l) as 
0.7 mag. As we see the red-shifted component of the ISO-Chal 217 outflow the outflow likely 
has a low inclination angle. Radial velocities of ~ -20 kms~^ and 40 kms~^ are measured 
for the red and blue lobes therefore tangential velocities of -55 kms~^ and 110 kms~^ are 
assumed for the purpose of this calculation. This is based on an outflow inclination angle of 
20°. 



Finally we also apply the method of Hartigan et al. ( 1995 ) to the blue-shifted outflow 



driven by ISO-Oph 102. The [SII]A6731 line is detected at a low S/N and its luminosity 
is estimated at 3.1 x 10^'' Lq. As the [SII]A6716 line is not detected there is no estimate 
of Ue for the outflow. However, as Uc for the [SII]A6731 line is closer to the value of Ug 
measured for CTT / BD outflows (than for the [OI]A6300 line). Equation 4 is less sensitive 
to the Uc / Ug dependence and we can make a reasonable estimate of the possible range of 
Ug for the ISO-Oph 102 outflow. Using a range of 10^ cm~^ to 10"^ cm~^ for Ug gives Mout at 
1.7-11.8 X 10"^'' MQyr^^. A value of Vtan = 78 kms~^ is assumed for this outflow based on 



the measured radial velocity of its FELs (Whelan et al. 2005) and on a outflow inclination 



of 30°. Phan-Bao et al. (2008) estimate the inclination angle of the ISO-Oph 102 disk at 63° 



to 66°. This suggests that the outflow inclination lies close to 30° as assumed. Rq is taken 



from Natta et al. (2004) where A„ is estimated at 3 mag. Note that Phan-Bao et al. (2008) 
measure Mout = 1-4 x 10~^ MqJt"-^ for the ISO-Oph 102 molecular outflow. It is reasonable 
that Mout for a molecular outflow be greater than the outflow rate in the underlying jet. 



Assuming that the jet is powering the molecular outflow (Downes & Cabrit 2007), the mass 



outflow rate of the molecular component will grow with time as the jet transfers increasing 
amounts of energy and momentum. 

Table 6. compares Mout and Mace for LS-RCrA 1, ISO-Chal 217 and ISO-Oph 102. As 
the typical width of a BD outflow is not known but reasonable estimates of outflow velocity 
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can be made from measured radial velocities, method B is more suitable at present as a tool 
for estimating Mout- An added advantage of method B is that Equation 4 is less sensitive to 
Ue which is useful for cases where Ug is not known and the [SII]A6731 line alone is detected. 
Method B has allowed us to limit Mout for all three objects to the range 10"^" to 10~^ 
M0yr~^. Quantities such as jet velocity and electron density must be better constrained in 
order to tighten this range. While the measured range for Mace is at present equally as large 
and is of order a few times 10^^° MQyr^^ up to of order 10~^, what can be concluded here 
is that Mout for the BDs in our sample is comparable to Mace- However, some important 
points to note are as follows. The number of objects in our sample is very small and the 
standard Mace versus Mout correlation for CTTSs shows a very wide scatter either side of the 
fitted least squares straight line ( Hartigan et al.|[r995 ). Importantly, there may also be some 
selection effects in our sample. As our observations are difficult and at the edge of what 
the VLT can do (the [OI]A6300 line emission is only just observable), one would therefore 
tend to pick up only objects with higher outflow rates (brighter [01] emission) all else being 
equal. In other words, in the plot (for CTTSs) mentioned above we might be selecting only 
the brighter (upper) points in the scatter. 

Difficulties related to measuring Maec in BDs should also be considered. The most 
common methods for deriving Mace are based on modelling of lines like Ha or Pa/? (Natta et 



al.|2004 ), on veiling measurements and on the flux of Call lines like the CaIIA8662 (Mohanty 



et al. 2005). However accretion rate estimates for BDs based on different lines (e.g. Halpha 



or Ca H&K) give very different results ( Mohanty et al.||2005 ). Individual lines (e.g. Ha) vary 
on time-scales of weeks ( Scholz et al.||2005 Nguyen et al.||2009 ) and thus the rates measured 
at a particular epoch may not reflect the long term accretion rate that is more appropriate 
for comparison with Mout- A second concieveable difficulty with these approaches is posed 



by any possible contribution to the Ha or Call lines from an outflow. For example Whelan 



et al. (2009) measure a signiflcant contribution to the LS-RCrA 1 outflow in the wings of 



the Ha line. Overall, further investigation is needed to conflrm if in BDs a larger fraction 
of infalling material is re-ejected in an outflow than in low mass stars. One theory of BD 
formation is that they originate from cores that never accrete enough material to enter a H 
burning phase. A mass outflow rate which is a signiflcant fraction of the infall rate could 
account for this. Future work on measuring Mout for young BD outflows should focus on 
increasing the number of objects known to have outflows, on obtaining better quality spectra 
and on imaging studies to allow outflow width and velocity to be derived. This work would 



also beneflt from expansion into the NIR regime (Podio et al. 2006). 



- 18 - 



Object 


Slit PA (°) 


He (cm ^) 








Te 


(10^ K) 




(10^ cm"^) 


LS-RCrA 1 


45 


1478 


(± 170) 


0.012 


(± 


0.005) 


35 


(± 1) 


1.2 


(± 0.7) 




63 


1300 


(± 140) 


0.004 


(± 


0.002) 


37 


(± 2) 


3.3 


(± 2.0) 




74 


1300 


(± 120) 


0.002 


(± 


0.001) 


41 


(± 1) 


6.7 


(± 2.3) 




82 


1500 


(± 200) 


0.002 


(± 


0.001) 


44 


(± 3) 


8.6 


(± 2.7) 


ISO-Chal 217 


B 


1900 


(± 700) 


















R 


5700 


(± 2300) 


















90 B 


600 ( 


± 200) 


















90 R 


2100 


(± 600) 

















Table 4: Physical conditions in the outflows of LS-RCrA 1 and ISO-Chal 217 as de- 
fined by ng, Xe, Te and iih- The electron density ng is estimated from the ratio of the 
[SII]A6716/[SII]A6731 lines. As explained above Xe, Tg and hh could only be measured in 
the case of LS-RCrA 1 and this was done using the BE technique. While estimates of Ug are 
reasonable when compared to measurements made for low mass protostellar outflows, values 
of Xg are far less and Tg, iih far greater than would be expected. Several factors make it 
currently very difficult to applying techniques such as the BE technique to BD spectra (refer 
to section 4.4) 



LS-RCrA 1 ISO-Chal 217 (red) ISO-Chal 217 (blue) ISO-Oph 102 



EW [OI]A6300 (A) 


-30.0 








-2 


EW [SII]A6731 (A) 


-6.4 




-0.7 


-0.3 


-0.2 


Ro 


19.2* 




16.59 


16.59 


14.57 


L [OI]A6300 (Lq) 


7.1 X 


10-^ 






3.1 X 10-^ 


L [SII]A6731 (Lq) 


1.4 X 


10-^ 


2.1 X 10~^ 


9.1 X 10-8 


3.1 X 10-^ 


Ug (cm~^) 


1300 




3900 


1250 


10^ -10^ 


in 


30 




20 


20 


30 


Vtan (kms"^) 


43 




55 


110 


78 


D (pc) 


130 




140 


140 


125 


Itan (") 


1.2 




1 


1 


1 



Hartigan et al. (1995). 



Table 5: Values used in the calculation of Mout from the method of 
As both lobes of the ISO-Chal 217 outflow are not resolved spectrally in the [OI]A6300 line 
it was not possible to measure their separate EWs so no values are given here. * Rq is the 



de- reddened magnitude and for LS-RCrA 1 is taken from Comeron et al. (2003). 
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Object 



Method 



LS-RCrA 1 



ISO-Oph 102 

ISO-Chal 217 Red flow 
ISO-Chal 217 Blue flow 



2.4 X 10^9 
6.1 X 10-1° 
2.0 X 10-1° 



A 

B [OI]A6300 
B [SII]A6731 



1.7-11.8 X 10-1° B [SII]A6731 

1.4 X 10-^ CO molecular outflowi 



3.1 X 10-1° 
1.8 X 10-1° 



B [S11]A6731 
B [S11]A6731 



^acc 



Method 



LS-RCrA 1 












2.8 X 


10- 


-10 


CaII(A8542)2 




10-1° 


-10" 


-9 


Optical Veiling'^ 




10 


-9 




CaII(A8662)4 




10" 


-10 




Ha 10% width^ 


ISO-Oph 102 












10 


-9 




Ha 10% width^ 




4.3 X 


10" 


-10 


J and K band spectra^ 


ISO-Chal 217 












1.0 X 


10' 


-10 


Ha emission^ 



Table 6: Measurements of Mout and Mace for the objects studied to date. For the BDs in our 



sample Mout is comparable to Mace- 1= 


Phan-Bao et al. 


(2008 


), 2= 


Comeron et al. 


(2003 


), 


3=Barrado y Navascues et al. (2004), 4= 


Mohanty et al. 


(2005 


, 5= 


Scholz & Jayawardhana 


(2006), 6=Natta et ah (2004), 7=Natta et ah (2006), 8= 


MuzeroUe et al. (2005) 
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5. Classical T Tauri-like Outflow Activity in the Brown Dwarf Mass Regime 



Much of what is now known about the launching of protostellar jets comes from inves- 
tigating the FEL regions of CTTSs. The first long-slit spectroscopic studies allowed their 



morphology, kinematics and position angle to be studied (Hirth et al. 1997) and progress in 



this area has been driven by high angular resolution techniques ( Ray et aL||2007 ). A notable 
feature of the FEL regions of CTTSs is that they often have two velocity components, a low 
velocity component (LVC) believed to trace a wide disk wind and a high velocity component 
(HVC) tracing the jet. Currently, the scale over which jets are collimated, the jet rotation, 
the jet excitation conditions and mass outfiow rate are all estimated from their FEL regions 
(Ray et al. 2007). While the distinctive FEL profiles associated with outfiow activity were 



first studied in the spectra of CTTSs (Hirth et al. 1997), Class I low mass protostars and 



Herbig Ae/Be stars are now also known to show strong FE. Studies have compared the FEL 
regions of the above objects to those of the CTTSs ( Boehm &: Catala|l994 Corcoran & Ray 
19981 [Davis et al.|[2003| [Garcia Lopez et al.|[2008| ). The earliest studies which identified FEL 



tracers of mass outflow in the spectra of BDs include Fernandez & Comeron (2001 ) and Bar 



rado y Navascues et al. (2004). Whelan et al. (2007) presented evidence which supported the 



hypothesis that BD outflows are scaled-down from CTT outflows. Here recent new evidence 
is summarised. As all the BDs investigated to date are optically visible it is assumed that 
they are the substellar mass counterparts to the CTTSs i.e. equivalent to Class II low mass 
sources. 

For the flve BDs whose FEL spectra have been investigated to date the [Ol]A6300 line is 
the dominant line. The [SII]A6731 line is detected in three cases (LS-RCrA 1, ISO-Chal 217 
and ISO-Oph 102) and only for LS-RCrA 1 is the [NI1]A6583 line detected at a reasonable 
S/N. Very faint [NII]A6583 emission is found in the spectrum of ISO-Chal 217. Interestingly 



for LS-RCrA 1 Whelan et al. (2009) detected both the blue and red-shifted outflow in the 



wings of the Ha line. This was the flrst time Ha was shown to originate in a BD outflow. 
As only blue-shifted FE is found this points to a dust hole in the LS-RCrA 1 disk, of an 
estimated size of ~ 7AU (refer to Whelan et al 2009 for further information). Spectro- 
astrometric studies of permitted emission lines have led to similar discoveries for the disk of 



CTTSs (Whelan et al. 2004 Takami et al. 2001). This is an important parallel between the 



disk-outflow systems of BDs and CTTSs. 

A second striking parallel between outflow activity in BDs and CTTSs is offered by the 
emergence of evidence of a LVC to BD FELs. Again this comes from the recent study of 
LS-RCrA 1. Although the LVC and HVC are not resolved spectrally for the LS-RCrA 1 



FELs, Whelan et al. (2009) argue that the large spread in radial velocity with the [OI]A6300 
line having the lowest velocity (-5 kms~^) and the [NII]A6583 line the highest (-22 kms~^ 
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with a wing extending to -75 kms~^) makes it clear that the [OI]A6300 hne is predominantly 
tracing a LVC component to the outflow and the [NII]A6583 line the HVC The [SII]A6731 
line shows evidence of both a LVC and a HVC. In the case of ISO-Oph 32 the low radial 
velocity of the [OI]A6300 emission line is also suggestive of the presence of a LVC to the 
outflow driven by this BD. An explanation for the low velocity based on a edge-on disk 
scenario is less likely given the lack of red-shifted emission (the 2MASS 1207-3932 FELs also 
have low radial velocities but both blue and red-shifted emission is seen) hence it is likely 
that the [OI]A6300 line has a strong contribution from a LVC. The conflrmation of a LVC 
to BD outflows would makes it very likely that many of the same processes active in the 
launching region of BDs and CTTSs and thus strongly links BD and protostellar outflows. 
However a large sample of young BDs must be investigated to obtain this conflrmation. 

As asymmetries are frequently observed in CTT outflows a third compelling piece of 
evidence supporting the continuation of CTT-like outflow activity into the BD regime, is 



the asymmetrical nature of the ISO-Chal 217 outflow. Hirth et al. (1994) first reported 



asymmetries in the jets driven by the CTTS RW Aur and DO Tau. For RW Aur the brighter 
red-shifted lobe has a radial velocity a factor of 2 slower than the blue-shifted lobe. In the 
case of DO Tau, while the blue-shifted lobe is the brightest it is slower than the red-shifted 



fiow by a factor of 2.3. In addition, the results of Hirth et al. (1994) reveal a difference in the 



electron density between both lobes of the flows. For RW Aur Ug is higher in the red lobe 



than the blue, however, for DO Tau the brighter lobe has the lowest electron density. Hirth 



et al. ( 1994 ) also list several other protostars for whom radial velocity asymmetries have been 



reported highlighting that this is a common feature of protostellar outflows. The discovery 
of such strong asymmetries so close to the driving source is important as it shows that these 
asymmetries are intrinsic to the central engine and suggests that their origin is in the jet 



launching mechanism itself. Ferreira et al. (2006) also discuss the origin of asymmetries in 



CTT jets. The asymmetries in the FELs of ISO-Chal 217 are remarkably similar to what 
has been well studied for RW Aur. In the spectra taken at both 0° and 90° the radial velocity 
of the red-shifted component is a factor of ~ 2 higher than that of the blue-shifted lobe, 
the red-shifted lobe is brighter and initial results suggest that the electron density of the 
red-shifted flow is also higher. Spectra taken along the actual outflow PA (estimated here at 
~ 200°) would allow a more detailed investigation of the ISO-Chal 217 outflow asymmetries. 

Finally we consider the existing evidence that BD outflows also have a molecular com- 
ponent. Although early observations of CTT jets concentrated on the optical regime ob- 



servations now extend into the near-infrared (Davis et al. 2003 Beck et al. 2008). Takami 



et al. ( 2004[ ) flrst reported a detection of a warm molecular component to the DC Tau jet 
(observed in the H2, v=l-0 S(l) line) and observations of H2 in the spectra of CTTSs are 
now not unusual. Note however that shocked H2 emission in outflows is much more common 
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where a source is in the embedded phase, i.e. Class - I ( Davis et al.||2001 ). While there has 
been very little investigation of BD outflow activity in the NIR, LS-RCrA 1 is known to show 
strong H2 emission at 2.122 /xm ( Fernandez fc Comer6n||2001 ). A recent significant discovery 
in the field of BD outflows is the observation of a CO J = 2 -1 molecular outflow driven by 
ISO-Oph 102 dPhan-Bao e'raL][2008| ). ISO-Oph 102 was the first BD shown by us to drive 
an outflow (Whelan et al. 2005). Blue-shifted FE was detected in its optical spectrum and 
found to be offset in a southerly direction. Phan-Bao et al. (2008) used the sub- millimeter 
array (SMA) to observe a small-scale molecular outflow driven by ISO-Oph 102. The total 
scale of the outflow is ~ 20" and it is aligned at a N-S PA. CTTSs are optically visible and 
therefore no longer embedded in their natal material. It is argued that molecular outflows 
originate as circumstellar material is entrained and driven outwards by an underlying jet 
(Downes & Cabrit 2007), therefore the most spectacular protostellar molecular outflows are 
driven by embedded protostars such as Class or Class I low mass stars (Arce et al. 2007). 
CO outflows are driven by the youngest CTTSs (Launhardt et al. 2009) and observations 
are rare. The existence of a molecular component to BD outflows is a fourth noteworthy 
similarity between BD and low mass protostellar outflows. 



6. Summary 

Currently our understanding of astrophysical outflows is focussed on observations of 
low mass protostellar objects. Although much can be gleaned from such observations a 
proper appreciation of the outflow phenomenon can only be reached through a complete 
investigation of how outflow activity compares for a range of driving sources, varying in 
mass and age. The new results included in this paper include the presentation of two new 
BD outflows, the unusual red-shifted asymmetry in the ISO-Chal 217 and the demonstration 
that Mout in BDs is comparable to Mace- Overall, the main aim of our BD project is to explore 
how substellar outflow activity compares to what is well documented for low mass protostars. 
As our sources are all optically visible our comparison is focussed on the CTTSs. A second 
advantage of this work is that it may add to the debate on the formation mechanism of BDs 
i.e. whether they form like low mass protostars or are ejected embryos from multiple systems. 
The later mechanism is thought to produce truncated disks and very little envelope. If BD 
outflows are long-lived this would argue against such a scenario although further investigation 
is required. 

Due to the faintness of the BD outflows observations are difficult and time consuming. 
However data gathered to date has led to some significant results. The most compelling 
pieces of evidence which support the continuation of CTT-like outflow activity into the BD 
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mass regime include, the detection of an outflow component to the LS-RCrA 1 Ha emission 
line and the suggestion of a possible disk dust hole, the identification of an LVC to the FELs 
of LS-RCrA 1 and ISO-Oph 32, the strong asymmetry in the ISO-Chal 217 outflow and the 
discovery of a molecular component to the optical outflows of LS-RCrA 1 and ISO-Oph 102. 

While these results are significant much further work is needed to provide a more com- 
plete picture of how BD and CTT outflows compare. This is especially true for any study of 
the physical conditions in BD outflows and BD mass outflow rates. A considerably increased 
sample is needed. This can be provided through further optical and NIR spectroscopic studies 
complemented by millimetre studies. Investigations must be extended to other wavelength 
regimes so that NIR atomic and molecular outflow tracers can be identifled and further CO 
outflows found. Finally optical/NIR images would allow the coUimation of the BD outfiows 
to be probed. Many of these projects are under-way or planned however due to faintness 
of the emission from a typical BD outfiows, it is hkely that long-slit /echelle spectroscopy 
combined with spcctro-astrometry will remain (in the short-term at least) the most effective 
way to investigate BD outfiow activity. 
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Fig. 3. — Spectro-astrometric analysis of the ISO-Cha217 [OI]A6300 line. The line has a 
blue and a red-shifted component. In the spectrum taken at 90° the blue-shifted component 
to the line appears as an extended wing. Offsets decrease at 90° and comparison with 
measurements at 0° constrain the outflow PA at ~ 208°. Offsets are a moving sum across 
the line and continuum position. The displacement in the line is measured in the continuum 
subtracted spectrum. The 1-a error in the line is marked by the dashed line. Finally as 
outlined in section 2 the line and continuum emission is summed in such a way that the 1-a 
error in both is comparable. 
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Fig. 4. — Spectro-astrometric analysis of the ISO-Cha217 [OI]A6363 line. Analogous to the 
[OI]A6300 hne, the hne is made of blue and red-shifted emission which spectro-astrometry 
shows to be offset in opposing directions. As expected the offsets are similar to what is 
measured for the [OI]A6300 line. Again comparison between measurements at 0° and 90° 
constrain the outflow PA at ~ 205°. As for all the ISO-Chal 217 FELs, offsets are a moving 
sum across the continuum and line position and the 1-a error in both measurements is 
comparable. 
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Fig. 5. — Spectro-astrometric analysis of the ISO-Chal 217 Ha line. Analysis was carried out 
in the same manner, as outlined for the [OI]AA6300,6363 lines. No signal from the outflow 
is detected in the Ha line. Emission is coincident with the source hence it is reasonable to 
assume that the bulk of the emission is coming from the accretion flow. For the Ha line the 
spectro-astrometric accuracy achieved at the line peak is much greater than in the wings or 
for the the continuum emission. Hence for all the Ha lines emission is summed so that the 
error in the line-wing and continuum is comparable. 
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Fig. 6. — Spcctro-astromctric analysis of the ISO-Chal 217 [SII]A6716 line. Again the line 
was observed at a PA of 0° and 90° degrees. The blue and red-shifted emission is more 
clearly separated spectrally than in the [OI]AA6300,6363 lines. In addition it is offset further 
from the BD position. Analysis of the [SII]A6716 line constrains the outflow PA at ~ 193°. 
This is smaller than but similar to the estimate from the [OI]AA6300,6363 hnes. Note the 
asymmetry in the red and blue lobes of the outflow in both velocity and intensity. 
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Fig. 7. — Spectro-astrometric analysis of the ISO-Chal 217 [SII]A6731 line. Analysis agrees 
with results for the [SII]A6716 line and constrains the outflow PA at ~ 200°. Again the red 
and blue lobes of the outflow are clearly separated and the asymmetry is especially clear. 
The difference between the intensity of the red and blue lobes as traced by the [SII]A6731 
line at 0° is striking. 
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Fig. 8. — Position velocity diagrams of the ISO-Chal 217 [SII]A6731 line region, comparing 
the line at slit PAs of 0° (left) and at 90° (right). For the PV diagrams the emission hne 
region was smoothed using an elliptical Gaussian filter as described above. Contours begin 
at 3 times the r.m.s noise and increase by factors of \/2- The black line delinates the BD 
position. The PV plots clearly support the spectro-astrometric results. The relative offset 
between the red and blue-shifted part of the line can be seen by eye in the plot for 0°. As 
this offset is a factor of 4 smaller at 90° it is not obvious in the corresponding PV diagram. 
A difference in brightness between the red-shifted components is also apparent with the red- 
shifted component at 90° appearing marginally brighter. As we are only detecting emission 
very close to ISO-Chal 217 and given that our slit width is 1 "we would not expect to see any 
significant variation in brightness with PA. The apparent difference here is due to a change 
in observing conditions between the acquisition of the two spectra with the seeing decreasing 
from 0'.'85 for the 0° spectrum to 0'.'7 at 90°. We measure the integrated flux of the red-shifted 
components to be comparable with the value being 0.2, 0.45 x 10^^^ crgs/s/cm^/A and 0.21, 
0.5 X 10~^^ ergs/s/cm^/A for the red/blue-shifted components at 0/90 degrees respectively. 
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Fig. 9. — [OI]A6300 line emission and related spectro-astrometric analysis for the young BD 
ISO-Oph 32. As reported for other BDs, the [OI]A6300 line was the only FEL found in the 
spectrum strong enough for spectro-astrometric analysis. Spectro-astrometry confirms the 
origin of the line in an outflow. The measured offset in the line is greater at 90° than at 0° 
and from this comparison the outflow PA is estimated at 240°. In the case of ISO-Oph 32 
the line and continuum regions were smoothed rather than summed. 
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Fig. 10.— Spectro-astrometry of the ISO-Oph 32 Ha line. As in the case of ISO-Chal 217 
offsets are a moving sum across the hne and continuum position and the hne and continuum 
are summed in such as way that the error in the hne-wings and continuum is comparable. 
The ± 1-a error is delinated by the dashed lines. As for all the sources discussed in this 
paper no signature of outflow activity is detected in the Ha line. 
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Fig. 11.— The Ha line and spectro-astrometric analysis for DENIS-P J160603.9-205644. 
DENIS-P J160603. 9-205644 is known to be actively accreting, however no evidence of outflow 
activity is found in our spectra. As is clear from the analysis shown here the Ho; emission 
is coincident with the source position and therefore must originate very close to the source. 
This analysis was carried out in the same manner as for all the other Ha lines. 
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